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Abstract. By employing the time-dependent Lanczos method, the nonequilibrium process of 
the half-filled one-dimensional extended Hubbard model under the irradiation of a transient laser 
>-^j ■ pulse is investigated. We show that in the spin-density-wave (SDW) phase, the antiferromagnetic 

C"| 1 spin correlations are impaired by the photoinduced charge carriers. Near the phase boundary 

Q , between the SDW and charge-density-wave (CDW) phases, a local enhancement of charge (spin) 

O ■ order that is absent in the original SDW (CDW) phase can be realized with proper laser 

frequency and strength. The possibility of restoration of spin orders from the CDW phase 
by optical means is discussed. 
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1. Introduction 

New insights into the dynamics of strongly correlated electron systems can be achieved by 
investigating their nonequilibrium processes, which demands the knowledge and understanding 
of the properties of excited states. The process have to be treated nonperturbatively if the 
system is driven far from the equilibrium. 

In experiments, one way to achieve nonequilibrium states is by applying various external 
electric fields, including pulses of electromagnetic field, which couple with charge degrees of 
freedom of the system. Well-known experimental facts on (quasi) one-dimensional (ID) materials 
include dielectric breakdown [TJ [2] , insulator to metal transitions and so on O HI [5] . 

In Mott insulators, we have to deal with the states above the Mott gap in order to induce 
charge carriers optically. It is well known that the ground states of Mott insulators have quasi- 
long-range antiferromagnetic spin orders. The charge carriers created by photons are, thus, 
expected to affect the spin correlations. It is desirable to obtain the knowledge of spin dynamics, 
as well as the charge's, during and after the doping. In many ID materials, it has been shown 
that besides the on-site Coulomb interaction U, nearest-neighbor (NN) interaction V has to be 
taken into account [UJ [7J- With increase of V, a phase transition from the Mott insulator to a 
charge-order phase can occur. The nonequilibrium dynamics near the phase boundary is also 
important. 



In this proceeding, we investigate the nonequilibrium process of a ID strongly correlated 
system under the irradiation of a laser pulse. The model we are working on is the half- filled 
extended Hubbard model, where the on-site, and NN interactions (U and V), are included. 
Its phase diagram in equilibrium is well-known and understood [HJ [TUl E] ■ Here we would 
like to reiterate some known properties of the model, concerning our work. In strong-coupling 
regime (large U(> 0)), with increase of V(> 0), a first-order phase transition from the spin- 
density-wave (SDW) to charge-density-wave (CDW) phase is realized when V reaches around 
[7/2. The transition can be understood as driven by the competition between the energetic costs 
for doublon generation and the benefits due to the attraction between doublons and holons. The 
SDW and CDW phases are recognized by algebraic decay of spin correlations, and a long-range 
(staggered) charge order, respectively. In both phases, the charge sector is always gapped except 
on the boundary. When spin degrees of freedom are concerned, gapless spin modes can be found 
in the SDW phase, while absent in the CDW phase. 

In the following sections, by studying time-dependent spin and charge correlation functions 
numerically, we show that spin correlations are suppressed by the photoinduced charge carriers in 
the SDW phase. Near the phase boundary from the SDW side, with proper laser frequency and 
strength, a sustainable charge order enhancement can be realized [12] but local spin correlations 
remains. Analogously, from the CDW side, the suppression of long-range charge order is 
accompanied with a local spin correlation enhancement. 

2. Model and Numerical Method 

The time-dependent Hamiltonian for the ID extended Hubbard model, where the external field 
is incorporated by means of the Peierls substitution for the hopping coefficients, can be written 
as 

H(t) = -t^(e iA Wci i(7 c m , ff + H.c.) 

+U £ (n i)t - 1) (n u -Vj+V^im-l) (n l+1 - 1) , (1) 

where cj a (cj )(J ) creates (annihilates) electrons with spin a at site i, n^ a = c| CT Cj j(T! rii = n^-f+ni^, 
th is the hopping constant. Particularly, with temporal gauge, we use a time-dependent vector 
potential A(t) to describe the laser pulse [T3] 

A(t) = A e- ( *-* o)2/2 ^ cos [oo pump (t - to)} , (2) 

where Aq controls the laser amplitude, which reaches its full strength at t = to; tj characterizes 
the duration time of light action. Note that, due to finite t^, the incoming photon frequency is 
broadened into a Gaussian- like distribution, with the variance of 1/^ around the central value 
w pump . Throughout this paper, we set th and t~^ 1 as energy and time units. 

Starting from the Schrodinger equation idip{t)/dt = H(t)ip(t), in order to obtain ip(t), we 
employ the time-dependent Lanczos method, which is originally described in Ref. [14], and 
followed by its applications in nonequilibrium dynamics of strongly correlated systems recently 
in Ref. [TS]. The basic idea is that we approximate the time evolution of \if}(t)) by a step- vise 
change of time t in small increments 5t. At each step, the Lanczos basis with dimension M is 
generated resulting in the time evolution 

M 

iw+st)) ~ e -^)«h/,(*)> ^^"ihmm)), (3) 
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(h) V=5.5 (A =0.60) 
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Figure 1. Time dependence of the spin-spin (left, from (a) to (d)) and charge-charge (right, 
from (e) to (h)) correlations as functions of distance (labeled by j) for 14-site lattice. The laser 
pulse with Gaussian magnitude modulation reaches its full strength at t = 12.5, as indicated by 
solid lines. Without exception, the pumping frequencies are set to match the resonance peaks 
of the optical absorption spectra. Parameters: St = 0.02, M = 100. In each subfigure, (a),(e): 
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where e\ and \(f>i), respectively, are eigenvalues and eigenvectors of the tridiagonal matrix 
generated in Lanczos iteration. As one of the advantages, the Lanczos method preserves the 
unitarity of the time-evolution operator at each time step, with relative low cost for a desired 
accuracy. The convergence check shows that for time length t around hundred, taking M = 30 
with 6t < 0.1 can provide adequate accuracy for the parameter regime we are working in. For 
even smaller 5t, M can be smaller. 

In the succeeding numerical calculations, periodic boundary conditions are employed to avoid 
possible boundary effects. We fix the parameters = 5, Iq = 12.5, and chose U = 10, where 
the first-order transition occurs around V ~ 5.1 [llj. 

3. Numerical Results 

In order to investigate the time evolutions of the spin and charge orders under the pules, we 
define the spin-spin and charge-charge correlations, respectively: 

S(j;t) = til J2(m\S*+j ■ SiW*)), (4a) 
i=Q 

C(j;t) = t^L J2(m\(ni +j - l)(n< - l)\m), (4b) 

i=0 

where L is the lattice size, and S, is the spin operator on site i. Intuitively, the antiferromagnetic 
spin and staggered charge orders are competitive. 

Figure [T] shows the results of correlations S(j; t), C(j; t), up to t = 50, for the system with 14 
lattice sites, where the largest distance between two sites is 7. Total crystal momentum P, as a 
good quantum number, is implemented. The size of the Hilbert space for P = is 841332. For 
2600 time step, it took around 30 hours for 12 cores. In succeeding discussions, we put more 
emphasis on spin. Detailed discussions on the charge correlations can be found in Ref. |12j . 

In Fig. [H four cases are shown: (a) and (e), deep into the SDW phase (V = 1); (b) and (f), 
the ground state is near the phase boundary, but set on the SDW side (V = 4.5); (c) and (g), 



(d) and (h), at the CDW side (V = 5.5) with different laser strength. Several observations can 
be made. First, notice that for the ground states (corresponding to t = 0) in the SDW phase, we 
can observe the gradually decay of spin correlations (Figs. [1(a), (b)) while no long-range charge 
order exists (Figs. H(e), (f)); in the CDW phase, the presence of a long-range charge order 
can be clearly recognized (Figs. H(g) or (h)), and the spin correlations largely vanish except 
for the NN one (See Figs. QJc) or (d), where S(j;0) < 1(T 4 for j > 1, and 5(1; 0) ~ 0.01). 
Second, for the case of V = 1.0, where the system is deep inside the SDW phase, introducing 
charge carriers into the system by the pulse cannot produce charge order (Fig. []Te)), while the 
spin correlations are severely suppressed except for the NN one (See Fig. QJa), where S(l; t) 
decreases from 0.14 at t = to 0.07 at t = 50, while some correlations with other distances 
even become negative). Notice that in this case, after the pulse, the total number of doublons 
(defined as iV^b = Si n i,t n «,i) ^ s around 2.2, corresponding to 30% doping [IB]. Third, we would 
like to pick up Fig. QJf), which shows that proper incoming laser pulse can enhance the charge 
order considerably in the SDW near the phase boundary, where we can observe that the NN 
correlation increases from 0.07 at t = to 0.22 at t = 50, together with a clear enhancement 
of other correlations in larger distances, which start from almost zero magnitudes at t = 0; 
meanwhile, the spin order does not suffer heavy suppression (See Fig. [ljb), where the NN 
correlation drops from 0.12 to 0.08). Analogously, from the CDW side close to the boundary, 
the pulse can suppress the charge order (Fig. QJg)) with local spin enhancement taking place 
at the same time (Fig. (He), where the NN correlation increases from 0.01 to 0.04). Finally, 
from Fig. [ljd) and (h), we show that strong enough laser pulse can destroy the charge order 
substantially without spin correlation enhancement. 

The results of the time-dependent correlations suggest intimate connections between spin 
and charge degrees of freedom in the photoexcited states. Interesting phenomena occur around 
the phase boundary, where local spin and charge orders can coexist. For a more quantitative 
analysis, we compare the spectra of systems with different values of V in the time-independent 
Hamiltonian, i.e., Eq. (pQ) but A(t) = 0. In order to describe the spin and charge orders which 
have finite spatial extensions, we define 

, L-l L c 

©sow = E E(- 1 )' S ^- • s * ( 5a ) 

c i=0 j=l 
-, L-l L c 

Ocdw = TF- E EC" 1 )* ^ - X ) fa - ^ ■ ( 5b ) 

c i=0 j=l 

Here L c is introduced as a cut off parameter for the correlation length, and we set L c = 5, 
estimated from the results of the density-matrix renormalization group (DMRG) on larger 
system size up to 30. The expectation values of the order parameter (Ocdw) and (Osdw) 
for eigenstates of a smaller 10-site system are shown in Fig. [21 

In Fig. [2] from (a) to (d), we can notice that for V = 1.0 to 5.5, states with relative large 
(Osdw) are always situated in the low-energy regime (say, E < 10). On the contrary, states with 
prominent (Ocdw) only appear there when the system is approaching the phase boundary. In 
Fig.[2]^g) for V = 4.5, around the resonance frequency of optical absorption spectrum (gjr ps 4.3), 
states with noticeable values of charge order can be found. This explains why a considerable 
charge order enhancement can be observed when the frequency of laser with proper strength is 
tuned to ojr (see Fig. Q^f), where ujr ~ 4.0 for 14-site lattice). However, in the complementary 
case, only the enhancement of the NN spin correlation is apparently visible, with rapid decrease 
with increase of distance. The same holds for L = 10, in spite of the fact that a state with 
excitation energy E ps 2.7 has a distinguished (Osdw) (shown by star in Fig. [2(d)). 
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Figure 2. The expectations of SDW (from (a) to (d)) and CDW (from (e) to (h)) order 
parameters of eigenstates for 10-site lattice with V = 1.0, 3.0, 4.5, and 5.5. The energy E is 
measured from the ground state. Only the data of states with P tot = <Stot = are shown, up to 
E = 60. Note that the largest values of the order parameters in (a) through (h) are marked by 
stars. 



In order to understand the different behavior of spin and charge correlations, and to further 
elaborate on the condition for the emergence of order enhancements induced by the laser pulse, 
we perform parameter-sweeping calculations in terms of uj pnmp and Aq for 10-site lattice. For 
a given pair of w pump and Aq, we carry out the time evolutions up to t = 52.5, then calculate 
the expectations of the SDW order , CDW order, and doublon number, noted as (OsDw)av, 
(0CDw)av) and (ndb)av) respectively, after averaging on the last 50 time steps (corresponding to 
time length At = 5). The results, presented in contour plots, are shown in Fig. [31 

In Fig. [3J we notice that for each V, there exist some frequencies with which the optical 
absorptions are most efficient, e.g., w pump « 7.5 in Fig. [3Ji) for V = 1.0, uj pmnp m 6.3 in Fig. [3J^j) 
for V = 3.0. These are resonance frequencies of optical absorption spectra ojr, which can be 
obtained from the imaginary parts of the dynamical current-current correlation functions. 

We can see that in order to excite the system effectively with weak fields, we need to tune 
the pulse frequency uj pump close to ojr. In general, in the SDW phase, accompanied with 
photodoping, spin correlations are suppressed (see the correspondence between the bright regions 
in Figs. [3](i) , (j), (k) for doublon generation, and the dark regions in Figs. [3(a), (b), (c) for 
(OsDw)av)- On the other hand, in the CDW phase, the destruction of charge order alone does 
not sufficiently bring back the spin order (see Fig. [3(1), compared with (d)). 

Things can be more subtle near the phase boundary. With Aq less than 0.1, a local 
enhancement of CDW order in the SDW phase can be induced accompanied with a suppression 
of spin orders, and vice versa in the CDW phase near the boundary. For example, see around 
Wpump ~ 4.3 in Figs. [3(g) and (c) for the case of V = 4.5, and w pump w 4.4 in Figs. [3^d) and 
(h) for the case of V = 5.5. Among them, the weakness of the SDW enhancement in Fig. [3]^d) 
can be attributed to the fact that the states with distinguished (Osdw) (around E w 2.7, see 
Fig. [2(d)) are in the off-resonance region, which makes it hard for them to be excited by the 
laser. 

Going to the high-field region, we observe that in Fig. [3(d), at w pump w 1.2, there is a well 
resolved stripe when Aq > 0.4, indicating an enhancement of the spin order. From the results 
of S(J; t) (not shown here), we notice that the spin correlations there extend themselves beyond 
NN. This may come from nonlinear multiphoton optical process |17j . Similar phenomenon for 
the case of charge order, can be found in Fig. [3(g), at uj pump « 2. Further investigation on the 
nonlinear effect is in progress. 
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Figure 3. Contour plots of final time-evolution results of SDW order (OsDw)av (first column), 
CDW order (C?cD\v)av (second column) and doublon number (last column) as functions of 
w pump and Aq, obtained by averaging on the last 50 time steps (time length At = 5), for 10-site 
lattice. w pump G [0.5,8]; Aq € [0.04,0.4] for the first three rows, and [0.08,0.8] for the last row. 
From top to bottom, V = 1.0, 3.0, 4.5, and 5.5, respectively. Here, we take St = 0.1, M = 30. 

4. Conclusions 

In summary, by studying time-dependent spin and charge correlation functions of the ID 
extended Hubbard model with applied laser pulse, we have shown that in the SDW phase, 
the spin order is suppressed due to the optical doping. Starting from the SDW side not far from 
the phase boundary, when the incoming laser frequency matches the optical absorption peak, 
a sustainable charge order enhancement can be realized. Similarly, starting from the CDW 
side, an enhancement of spin correlations, though more localized, has been observed. In the 
off-resonance region, we find more extended recovery of spin correlations which may come from 
nonlinear effects. 
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